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<N . Abstract 

We demonstrate that strong suppression of the relative production rate (d+Au) / (p+ 
p) of inclusive high-p^ hadrons at forward rapidities observed at RHIC is due to 
parton multiple rescatterings in nuclear matter. The light-cone dipole approach- 
^ ■ based calculations are in a good agreement with BRAHMS and STAR data. They 

also indicate a significant nuclear suppression at midrapidities with a weak onset 
of the coherence effects. This prediction is supported by the preliminary d + Au 
data from the PHENIX Collaboration. Moreover, since similar suppression pattern 
Qh| is also expected to show up at lower energies where effects of parton saturation are 

not expected, we are able to exclude from the interpretation of observed phenomena 
models based on the Color Glass Condensate. 
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1 Introduction 

rn ; 

t— ( ■ Spectra of high-p^ hadrons produced in nuclear collisions at large forward rapidities are 
promising tool to study partonic degrees of freedom in nuclei. Strong nuclear suppression of 
the spectra observed by the BRAHMS [H [2] and STAR [3] Collaborations in deuteron-gold 
j> \ collisions at the Relativistic Heavy Ion Collider (RHIC) was a tempting invitation for the 
parton saturation [U [5] or the Color Glass Condensate (CGC) j6] motivated phenomenology 
[7] as its most natural explanation. 

According to these models the parton coherence phenomena may reveal itself already at 
RHIC energies showing up first in the wave function of heavy nuclei. Kinematically most 
favorable region to access the strongest coherence effects is the fragmentation region of the 
light nucleus 1 colliding with the heavy one 2. At large x\ (large Feynman at forward 
rapidities) one can reach the smallest values of the light-front momentum fraction variable 
x 2 = X\ — xp (2 x 10~ 4 <x 2 < 10~ 3 in the RHIC kinematic range). 

Quite unexpectedly, the same nuclear effects occur not only at forward rapidities [H [21 [3] 
but also in the large pt region at midrapidity [8] where effects of coherence are not important. 
The covered interval of x 2 ^0.01 goes too far beyond the region where the CGC is valid. 

In [HI HO] it was shown that for any large-Xi reaction considerable nuclear suppression 
comes from the energy conservation at the level of projectile partons undergoing multiple 
rescatterings in nuclear medium. It was also demonstrated [9] that large-^i suppression is a 
leading twist effect, violating QCD factorization, a basic ingredient of the CGC-based models. 

Analysis of nuclear suppression based on the multiple parton rescatterings leads also to a 
new type of scaling: the same nuclear effect are expected at different energies and rapidities 
corresponding to the same value of x\ (xf at forward rapidities) [9|,[T0]. The most straight 
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forward prediction of the Xi-scaling is that similar nuclear effects must also show up at lower 
cm. energy y/s. Here the onset of coherence effects is much weaker and so there is much less 
room for explanation of strong nuclear suppression in terms of the CGC. 

Another consequence of this scaling is that in the RHIC energy range similar nuclear 
effects must also show up at midrapidities provided that the corresponding values of px of 
produced hadrons reach the same value of X\ as at forward rapidities. This prediction is 
confirmed by the preliminary data on neutral pion production in d + Au collisions measured 
recently by the PHENIX experiment |8J showing an evidence for the nuclear suppression at 
rather large px > 8 GeV. This and new 2008 d + Au high-statistic data may provide another 
test of our approach. 

2 High-/>r hadron production: Sudakov suppression, pro- 
duction cross section 

Let us recall that in the limit X\ — *■ 1 (xp — > 1 at forward rapidities) gluon radiation in 
any pQCD-driven hard scattering is forbidden by the energy conservation. For uncorrelated 
Poisson distribution of radiated gluons, the Sudakov suppression factor, i.e. the probability 
to have a rapidity gap Ay = — ln(l — x\) between leading parton and rest of the system, has 
a very simple form: S{x\) = 1 — x\ [9]. 

Suppression at X\ — > 1 can thus be formulated as a survival probability of the large 
rapidity gap (LRG) process in multiple interactions of projectile valence quarks with the 
nucleus. Every additional inelastic interaction of the quarks contributes an extra suppression 
factor S(xi). The probability of an n-fold inelastic collision is related to the Glauber model 
coefficients via the Abramovsky-Gribov-Kancheli (AGK) cutting rules |TT] . Correspondingly, 
the survival probability at impact parameter b reads 
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where T^(b) is the nuclear thickness function. 

At large px, the cross section of hadron production in d + A (p + p) collisions is given 
by a convolution of the distribution function for the projectile valence quark with the quark 
scattering cross section and the fragmentation function 
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parametrization from [T2]. Fragmentation functions were taken from 

As first shown in [9j [10] the effective projectile quark distribution correlates with the 
target. So interaction with the nuclear target does not obey the factorization. Main source of 
suppression at large px comes from multiple soft rescatterings of the quark in nuclear matter. 
Summed over multiple interactions, the quark distribution in the nucleus reads 



/^k^iiSr) = C f q/N {xi,q T ) 



Jd 2 b 


e —xicr e ffTA(b) 




(1- 


Xi) j d 2 b 


1 - 







(3) 



where effective cross section a e ff = <7 e //(pT, s) = t^tt^ has been evaluated in The 
normalization factor C in Eq. ([3]) is fixed by the Gottfried sum rule. 

The cross section of quark scattering on the target da[qA{p)]/ dpqxdri in Eq. (jSJ) is calculated 
in the light-cone dipole approach [HI [15]. We separate contributions with different initial 
transverse momenta and sum over different mechanisms of high-p^ hadron production. Details 
can be found in [H]. 

Let us note that in the RHIC energy range and at midrapidity correct description of 
hadrons with small and moderate pt can be achieved only if the above calculations incorpo- 
rate production and fragmentation of gluons[T6]. Consequently, the cross section for hadron 
production, Eq. (j2J), should be extended also for gluons with corresponding distribution func- 
tion, parton scattering cross section and the fragmentation function. Including multiple par- 
ton interactions, the gluon distribution in the nucleus is given by the same formula as for 
quarks (see Eq. (j3J)), except cr e ff, which should be multiplied by the Casimir factor 9/4. 



3 Comparison with data 

In 2004 the BRAHMS Collaboration [TJ observed a significant nuclear suppression of negative 
hadrons produced at r\ — 3.2. Much stronger onset of nuclear effects was observed later on by 
the STAR Collaboration [3] for ir° production at pseudorapidity r\ = 4.0. Both measurements 
are plotted in the left panel of Fig. [U A huge difference in nuclear suppresion factor at 
different r\ is due to the energy conservation and reflects much smaller survival probability of 
the LRG in multiple parton interactions at larger xi [U [TO] . 
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Figure 1: (Left) Ratio of negative hadron and neutral pion production rates in d + Au and p + p 
collisions as function of px at r] = 3.2 and r\ = 4.0. Data are from the BRAHMS pQ and STAR 
Collaborations [3], respectively. (Right) Model predictions for the ratio Rd+Au(pr) for production 
of 7T° mesons at y/s = 200 GeV and different values of rj changing from 3 to 4. 



To demonstrate different onsets of nuclear effects with increasing pseudorapidity we present 
in the right panel of Fig. [TJ our calculations for the nuclear suppression factor at different fixed 
values of 77. Changing the value of r\ from 3.0 to 4.0 leads to a rise of Rd+Auipr) by a factor 
of 2 M- 
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Figure 2: (Left) Ratio, Rd+AuiPr) for identified particles produced in d + Au and p + p collisions 
at rj = 3. The data are from the BRAHMS Collaboration [2\. (Right) Predictions for the ratio 
Rd+AuiPr) for 7T° at different rj and having the same exp^/y^s. 



The BRAHMS Collaboration has recently reported a new measurements [2] on production 
of positively charged pions and kaons at 77 = 3.0 in d + Au collisions confirming suppression 
pattern they found in 2004 for the negative particlespQ. Their recent data are plotted on the 
left panel of Fig. [2] together with our model predictions. 

The calculations of Rd+Auipr) of neutral pions at yfs = 200, 130 and 62.4 GeV shown on 
the right panel of Fig. [2] reveal approximate xi(xF)-scaling at RHIC energy range, i.e. the 
same nuclear effects at values of i] and yfs corresponding to the same value of x\. 

Generalization of the ^-scaling from the forward region to midrapidity is studied on Fig. [31 
The only difference to the previous analysis is that the same value of X\ at midrapidity as 
that in the forward region requires substantially higher hadron transverse momenta. On 
the left panel of Fig. [3] our predictions for the nuclear suppression factor of 7r° produced 
in d + Au collisions at midrapidities are confronted with the recent data of the PHENIX 
Collaboration [8J . 

Here the thin dashed line corresponds to the case when multiple parton rescatterings are 
not taken into account. The calculations with inclusion of multiple parton rescatterings are 
presented by the thin solid line. At moderate px £ (3, 7) GeV the model underestimates the 
data. However, quite a strong onset of nuclear suppression at large px is not in a disagreement 
with corresponding experimental points. At px = 25 GeV we expect Rd+Auipr) ~ 0.9. 

Due to the transition between the regimes with (small px) and without (large px) onset of 
coherence effects in the RHIC energy range calculations at 77 = are very complicated. One 
can deal with this situation relying on the light-cone Green function formalism [TTl [TEJ [19] 
but the integrations involved become too complicated. To simplify the situation we have used 
instead corrections for finite coherence length. Following the procedure described in [16] we 
have used linear interpolation performed by the means of so-called nuclear longitudinal form 
factor. Such a situation is shown by the thick solid and dashed lines on Fig. [3] corresponding 
to the case with and without inclusion of the multiple parton rescatterings, respectively. One 
can see that this correction brings the model predictions to a better agreement with the data 
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Figure 3: (Left) Ratio Rd+Auipr) as a function of pr for production of 7r° mesons at \fs = 200 GeV 
and rj = vs. data from the PHENIX Collaboration [8]. Thin solid and dashed lines represent the 
predictions calculated in the limit of long coherence length. Thick solid and dashed lines include 
corrections for the finite coherence length. (Right) The same as Fig. in the left panel but for the 
ratio R p+ Au{pt)- 

at moderate pr- 

On the right panel of Fig. [3] we also present model predictions for the ratio R p+ a u as a 
function of pt- Compared to d + Au system study of nuclear effects in p + Au minimizes the 
isospin effects. At px = 25 GeV we predict R p +a u ~ 0.93. 

4 Summary and conclusions 

In the present paper we have analyzed consequences of the X\ (x^)-scaling of the nuclear 
suppression factor R p (d)+Au of high-p T hadrons at RHIC. 
The new results are: 

• According to the Xi-scaling, considerable nuclear suppression at large x\ is expected for 
different kinematic regions : 

- production of high-p^ hadrons at forward rapidities. 

- production of high-pr hadrons at smaller rapidities and smaller energies. 

- productions of hadrons with very large at midrapidity. 

• Using simple formula ([3]) adopted from [9] and based on the Glauber multiple interaction 
theory and the AGK cutting rules, we have calculated high-p^ hadron production at 
midrapidity and found quite a strong nuclear suppression. This observation does not 
contradict to the recent measurements of the PHENIX Collaboration [8]. 

• In order to avoid the isospin effects, we have also studied large-pr neutral pion produc- 
tion in p + Au collisions. With the same input, we predict (see the right panel of Fig. [3]) 
for the first time quite a strong nuclear suppression, R p+ a u = 0.93 at px = 25 GeV. 

• In the RHIC kinematic region, investigation of large- a; hadron production in p(d) + Au 
collisions at midrapidities represents the baseline for verification of different phenomeno- 



logical models, important. At high-p T the data cover region of x 2 ~ 0.05 — 0.1 where 
effects of coherence are negligible allowing to exclude the CGC-based models from in- 
terpretation of observed nuclear suppression. 
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